Entanglement of strongly interacting low-dimensional fermions in metallic, superfluid 

and antiferromagnetic insulating systems 



o 
o 

(N 

o 
O 

(N 



u 



X3 
O 

(N 
> 

O 

(N 

o 
o 



o 



X 



V. V. Franga and K. CapelleQ 

Departamento de Fisica e Informdtica, Instituto de Fisica de Sao Carlos, 
Universidade de Sao Paulo, Caixa Postal 369, 13560-970 Sao Carlos, SP, Brazil 

(Dated: February 6, 2008) 

We calculate the entanglement entropy of strongly correlated low-dimensional fermions in metallic, 
superfluid and antiferromagnetic insulating phases. The entanglement entropy reflects the degrees 
of freedom available in each phase for storing and processing information, but is found not to be 
a state function in the thermodynamic sense. The role of critical points, smooth crossovers and 
Hilbert space restrictions in shaping the dependence of the entanglement entropy on the system 
parameters is illustrated for metallic, insulating and superfluid systems. The dependence of the 
spin susceptibility on entanglement in antiferromagnetic insulators is obtained quantitatively. The 
opening of spin gaps in antiferromagnetic insulators is associated with enhanced entanglement near 
quantum critical points. 
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I. INTRODUCTION 

Entanglement, one of the most surprising predictions 
of quantum theory, has recently received much atten- 
tion in the context of quantum information theory and 
quantum computation. In quantum information theory, 
entanglement is identified with nonseparable wave func- 
tions, and referred to as "quantum correlation", as op- 
posed to " classical correlation" . In quantum many-body 
physics and quantum chemistry, the expression "corre- 
lation" is used in a much more restricted sense, mean- 
ing those consequences of the particle-particle interac- 
tion that arise beyond the mean-field approximation. 
Both meanings of correlation are conceptually distinct: A 
Hartree-Fock wave function (a single determinant), e.g., 
is entangled in the quantum-information sense, but not 
correlated in the many-body sense. Even a wave function 
for perfectly noninteracting and spatially widely sepa- 
rated particles can be entangled (as illustrated by the 
Einstein-Podolsky- Rosen paradox), whereas correlations 
in the many-body sense vanish for strictly noninteracting 
particles. 

The present paper investigates correlations in the 
quantum information sense, in a model in which correla- 
tions in the many-body sense are known to be strong: the 
fermionic Hubbard model. The one-dimensional Hub- 
bard model in an external magnetic field has the Hamil- 
tonian 

H = ^(cLcj+i,,, + H.c.) + UJ2 4\(^^l4i'^il 

i.a i 

-f^J2''l'^''"' ^ \ ^(4^'^i^ - clic^i), (1) 

where t describes hopping between neighbouring sites, U 
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is the on-site particle-particle interaction, h is a mag- 
netic field, /i the chemical potential, and the operators 
Cio- satisfy fermionic commutation relations |jy]. This 
Hamiltonian has a long tradition as archetypical model of 
strongly interacting electrons in solids comprising 
metallic and insulating phases at repulsive interactions, 
and a phase with superconducting correlations at attrac- 
tive interactions. More recently it has also been used 
to describe fermionic atoms confined in optical lattices 
^ IE 111' In the absence of spatial inhomogeneity this 
model has an exact solution in terms of the Bethe-Ansatz 
(BA) 0) 0> ; which reduces the problem of finding the 
many-body ground state and its energy to solving a set 
of numerically tractable coupled integral equations. 

Differently from the Heisenberg model, which is more 
widely used in investigations of entanglement, the Hub- 
bard model also accounts for the effect of itineracy of 
particles and possible superconductivity. What are the 
effects of such complications on entanglement measures 
and entanglement witnesses? 

Entanglement measures for the Hubbard model were 
proposed and calculated, e.g., in Refs. 1^ IT(1 ITU 
il^ il4j . Specifically, our work builds on, and extends, 
that of Refs. |3 and on the local von Neumann en- 
tropy Tij , or single-site entanglement £ 1, 0, 0| of the 
Hubbard model. Gu et al., in Ref. 8], numerically study 
entanglement on finite Hubbard clusters. For L = 70 
sites, they obtain the entanglement entropy as a function 
of U precisely at haf filling {n = 1), whereas for smaller 
clusters of up to L = 10 sites they also obtain results for 
n 7^ 1. Larsson and Johannesson Q consider the oppo- 
site limit, L — > cxD, in which they obtain the entanglement 
entropy of noninteracting particles in a half-filled band 
{n = 1,U = 0) and of particles with strongly attractive 
interaction in a half-filled band (n = l,U ^ — oo) as a 
function of magnetic field h, and of particles with repul- 
sive interaction in zero magnetic field (h — 0,U > 0) as a 
function of chemical potential /i |l(i| . The main interest 
of Refs. 1^ and was in the behaviour right at critical 
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points, where the entanglement or its derivatives were 
found to be strongly enhanced. 

In the present paper, we use an efficient evalua- 
tion of the Bethe-Ansatz integral equations, recently 
developed in the context of density-functional theory 
HIEollEIH) in order to explore the entire n{^)-U-h 
phase diagram, both at and away from critical points 
|20|. Special cases, known from the earlier work described 
above, are recovered, and complemented by new infor- 
mation on the behaviour in other regions of the phase 
diagram. We also evaluate an observable not considered 
in Refs. Q and 0, namely the spin susceptibility, known 
to be an entanglement witness pH l23j . Our results 
on the spin susceptibility have a direct bearing on the 
recent reanalysis 24] of earlier experiments |2?t| of the 
entanglement in antiferromagnetic cuprate ladders, and 
also provide a new perspective on the enigmatic spin gap 
observed in the normal state of cuprate superconductors. 



II. ENTANGLEMENT ENTROPY ACROSS THE 
PHASE DIAGRAM 

Starting point of the present analysis is the expression 
of £ in terms of the particle density (or filling factor) 
n = n-^ + ni = (c||Ci|) + (cl^Cii), the magnetization (or 
spin density) m = (ri| — ri|)/2, and the particle-particle 
interaction U, 
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where e = E{n,m,U) / L is the ground-state energy per 
site. 

Equations equivalent to JSJ can be found in Refs. 0,01 j 
where they are evaluated by approximating the energy 
derivative dE{n,m,U) / dU in certain special limits 9], 
or for systems with a few lattice sites in order to 
study £{n,m,dE/dU) near critical points. Our interest 
here is in the behaviour of the entanglement entropy also 
away from critical points and special limits. To this end 
we must have access to the full function £{n, m, dE/dU). 
Below, we obtain this function from numerical solution 
of the BA equations for e{n,m,U). In the context of 
applications of density-functional theory to spatially in- 
homogeneous Hubbard models y,aLJ,UJ,LJ we have 
recently obtained such numerical solutions on a dense 
mesh of values of n, m and U. 

Figures n and m illustrate the resulting dependence of 
the entanglement entropy on the particle-particle inter- 
action and the particle density. Increase of \U\, which 
makes the wave function more correlated in the many- 
body sense, decreases making it less entangled in the 



quantum information sense. Conversely, for noninteract- 
ing particles {U = 0) the entanglement entropy reaches a 
maximum. This behaviour of £{U) is due to a nontrivial 
restriction on the Hilbert space of the one-band Hubbard 
model, which permits at most one particle of each spin 
on each site. This restriction remains effective even for 
U = (no correlations in the many-body sense) and im- 
plies strong correlations (entanglement) in the quantum 
information sense. We expect that similar Hilbert-space 
restrictions (e.g., due to conservation laws or topology) in 
more complex systems lead to other examples of counter- 
intuitive relations between both notions of "correlation" . 

The qualitative change at the critical points U = and 
n — 1 was predicted in Refs. but Figs. Q and 13 re- 

veal that even away from critical points the dependence 
of £{n,m,U) displays considerable structure. For posi- 
tive U (repulsive interactions) the entanglement curve as 
a function of U is almost flat, except at the critical point 
n = 1, where a pronounced U dependence is observed, 
and the curve drops from its theoretical maximum £ — 2 
to below the values observed for n = 0.5. This differ- 
ent behaviour at n = 1 reflects the Mott metal-insulator 
transition. The rapid drop of the entanglement entropy 
in the insulating phase is due to the freezing of electronic 
degrees of freedom at n = 1 and large positive U, asso- 
ciated with the transition from the physics of itinerant 
electrons to that of localized antiferromagnetically cou- 
pled spins. At n 7^ 1 the electrons remain itinerant and 
the translational degrees of freedom are not frozen out. 

We observe in Figs. n]and|llthat curves corresponding 
to different values of the system parameters cross, show- 
ing that physically distinct states can give rise to the 
same value of the entanglement entropy, which is thus 
not a faithfuU state function reflecting the microscopic 
structure of the ground state. 

At negative U (attractive interactions) a rapid drop in 
£{U) is observed for all fillings. This drop is again due to 
freezing of degrees of freedom, but at negative U the rel- 
evant physics is the crossover from weakly coupled pairs 
in a BCS-like state to strongly coupled localized dimers 
in the Bose-Einstein (BE) limit. The number of degrees 
of freedom is reduced by a factor of 2 in the formation 
of tightly bound pairs. Accordingly, the numerically de- 
termined ratio £{U = 0)/£{U -oo) « 2. The en- 
tanglement entropy thus serves as a marker not only for 
quantum criticality |^ ^2d\ but also for smoother 

crossovers. The strong dependence of entanglement on 
the size of the pairing interaction must be taken into ac- 
count in investigations of entanglement and qubits in the 
superconducting state 27, 28, 29, 30J. 

From the point of view of density-functional theory, in 
particular the Hohenberg-Kohn theorem, the fundamen- 
tal variables, determining the value of any observable, are 
the densities n and m, and the interaction U. From the 
point of view of experimental manipulation of entangled 
particles in a device, it is, however, much more conve- 
nient to cast the results in terms of variables that are 
more directly controlled in the laboratory, such as the 
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FIG. 1: Local entanglement entropy for /i = as a function 
of particle-particle interactions, for different particle densities 
n. 
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FIG. 2: Local entanglement entropy for ft = as a function 
of particle density, for different interactions U . 



applied magnetic field h or the chemical potential /i. 

Given the function £{n, m, dE/dU), the relations fi — 
dE{n, m, U)/dn and h = dE{n, m, U)/dm can be used to 
numerically construct the dual function £{fj,,h,U). Nu- 
merical results for £{^{n), h — 0,U) and ii{n) are plotted 
in Fig. O which shows that the entanglement entropy is 
a very rapidly varying function of the chemical potential. 
The behaviour at [/ < is similar to that at [/ = 0, 
but for U > a gap opens at n = 1 in the energy spec- 
trum, the chemical potential becomes ill-defined, and the 
entanglement entropy displays a pronounced derivative 
discontinuity. 

Figure 01 illustrates the magnetic-field dependence of 
the entanglement entropy. For all the curves termi- 
nate at the magnetic field corresponding to magnetic sat- 
uration. For U > 0, the entanglement entropy is maxi- 
mal at /i = and decreases monotonically with increasing 
magnetic field, refiecting increased magnetic order. For 
U < 0, the entropy only starts to depend on \h\ for fields 



FIG. 3: Local entanglement entropy at /i = as a function 
of chemical potential (upper panel), and chemical potential 
as a function of density (lower panel), for different values of 
the interaction U. Note that both panels share a common fj, 
axis, permiting one to reconstruct the function £{n). 




FIG. 4: Local entanglement entropy at n 
of magnetic field h for several values of U. 
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above the spin gap, and displays a round maximum at 
some \h\ > 0. 



III. SPIN SUSCEPTIBILITY 

As seen above, the behaviour of the entanglement pro- 
vides direct information about the degrees of freedom 
available for storing and processing information in sys- 
tems of strongly interacting fermions. Other observ- 
ables are also strongly affected by changes in the en- 
tanglement. We here focus on the spin susceptibility, 
(i) because this quantit y is an example of an entan- 
glement witness 0, , and (ii) because of recent 
claims p3| that anomalies earlier observed in the 
spin susceptibility of the antiferromagnetic (AFM) in- 
sulator Cu(N02)22.5H20 are due to enhanced entangle- 
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ment. This claim is based on the observation that a puz- 
zhng and hitherto unexplained (251] drop in the spin sus- 
ceptibility begins to develop at approximately the same 
temperature where entanglement is predicted to increase 

However, as already indicated in the experimental pa- 
per p^l, many different explanations for the observed 
anomaly are conceivable. Moreover, in recent analysis of 
experiments on LiH0xYi_xF4j with predominantly ferro- 
magnetic interactions, precisely the opposite behaviour, 
namely an increase in the spin susceptibility as entangle- 
ment laecomes stronger, was argued to occur jSlj . 

It is thus important to verify the connection between 
the susceptibility and entanglement in a model stripped 
off all unessential complications, but still complex enough 
to encompass both antiferromagnetism and entan- 
glement. The spins in Cu(N02)22.5H20 are arranged 
in chains, so that a description in terms of the one- 
dimensional Hubbard model, which at n = 1 and U —* oo 
is equivalent to the AFM Heisenberg model with J ~ 
4t^/U, becomes possible. 

The zero-temperature spin susceptibility of the Hub- 
bard model can be related to the ground-state energy 
via = d^E{n,m, U) / dm?\m=o-, and may be obtained 
numerically from the Bethe Ansatz equations j32| as a 
function of U. From Fig. ^we see that for [/ > and 
n — 1 the function S{U) is invertible. Hence, we can 
numerically construct the function x~^{U{£)) — x~^{£)- 
For n = 1 (consistent with the fact that the experiments 
were done on insulators), this function is plotted in Fig. 
which clearly shows that the susceptibility decreases 
with increasing entanglement entropy. Hence, in a min- 
imal (but still realistic) model of antiferromagnetic in- 
sulators the susceptibility indeed drops sharply as the 
system's wave function becomes more entangled. This 
result strongly suggests that the observation of a drop in 
the spin susceptibility in Cu(N02)22.5H20 ]^ as entan- 
glement grows (2i| is not a coincidence, or an artifact. 

These Hubbard-model calculations, and according to 
also the experiments of Ref. show that en- 

hanced entanglement in AFM systems produces a sharp 
drop in the spin susceptibility. The Hubbard-model cal- 
culations reported in the first part of this Letter, as well 
as related work in |^ |^, ^| , show that in the proximity 
of quantum-critical points the entanglement entropy or 
its derivatives can be very strongly enhanced. By putting 
these two pieces of information together we infer that in a 
strongly entangled state the spin susceptibility of antifer- 
romagnetic systems is suppressed near quantum critical 
points. 

Interestingly, there is another important class of ma- 
terials that display this type of behaviour, although it is 
not normally discussed in these terms: The parent com- 
pounds of high-temperature superconductors are AFM 
cuprates, for which much experimental evidence points 
to the existence of a quantum-critical point. The spin 
susceptibility of these systems displays a puzzling drop 
at low temperatures (see, e.g.. Fig. 4 and p. 16004 of 
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FIG. 5: Spin susceptibility at n = 1 and /i = as a function 
of the entanglement entropy. 



1331 ■ and the recent re-analysis of those data in Fig. 4 
of |34|). This drop is very similar to the one observed 
in the AFM spin chain Cu(N02)22.5H20, and suggests 
that the opening of a spin gap in high-temperature su- 
perconductors is also be associated with entanglement, 
persisting at relatively high temperatures. 



IV. CONCLUSIONS 

In summan^ we have confirmed the previous obser- 
vation 1^ 0, 0, that the entanglement entropy, 
or its derivatives, are strongly enhanced near quantum 
critical points. Our Bethe- Ansatz-based treatment ad- 
ditionally reveals that (i) similar behaviour is observed 
at crossovers, (ii) even away from critical points or 
crossovers, each physically different phase leaves its dis- 
tinctive mark on the entanglement entropy, (iii) but in 
spite of this, the entanglement entropy is not a state 
function in the thermodynamic sense. It does, however, 
(iv), provide detailed information about which degrees 
of freedom are available for storing and processing infor- 
mation in each of the possible phases of the system, (v) 
Hilbert-space restrictions can lead to counterintuitive re- 
lations between the concept of correlation, as employed in 
many-body physics, and that of entanglement, employed 
in quantum-information theory. 

In addition to the entropy, we also extract another ob- 
servable closely related to entanglement: the spin sus- 
ceptibility. We find numerically that in AFM systems 
enhanced entanglement sharply reduces the s pin suscep- 
tibility, in agreement with recent reanalysis 24] of ear- 
lier experiments 25j but in sharp contrast to what was 
found in ferromagnetic systems [IJ. This observation 
suggests that the enigmatic spin gap, observed above 
the critical temperature in cuprate superconductors, may 
also be a consequence of enhanced entanglement due to 
proximity of a quantum critical point. Note that high- 
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temperature entanglement in cuprates was suggested on 
different grounds also in |33|. 
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